Spin Accumulation in the Extrinsic Spin Hall Effect 
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The drift-difFusion formalism for spin-polarized carrier transport in semiconductors is generalized 
to include spin-orbit coupling. The theory is applied to treat the extrinsic spin Hall effect using 
realistic boundary conditions. It is shown that carrier and spin diffusion lengths are modified by 
the presence of spin-orbit coupling and that spin accumulation due to the extrinsic spin Hall effect 
is strongly and qualitatively influenced by boundary conditions. Analytical formulas for the spin- 
dependent carrier recombination rates and inhomogeneous spin densities and currents are presented. 
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In the presence of spin-orbit coupling, either due to 
impurities or due to host lattice ions, carriers of op- 
posite spins tend to scatter into opposite directions. 
With an electric field induced (longitudinal) motion un- 
der bias, the spin-orbit scattering results in a transverse 
spin current and spin accumulation, as first predicted by 
D'yakonov and Perel' and later revisited by oth- 

ers . This effect .which is now called the extrinsic spin 
Hall effect (SHE) 4], has been recently demonstrated ex- 
perimentally in n-GaAs and n-InGaAs thin films ^ and 
in two-dimensional electron gas confined within (110) Al- 
GaAs quantum wells ^ . The signature of the effect is op- 
posite spin accumulation at the edges of the sample, with 
spin polarization perpendicular to the transport plane. 

This paper has two goals. First, we present a formal- 
ism for carrier drift and diffusion in inhomogeneous spin- 
polarized semiconductors in the presence of spin-orbit 
coupling and spin-dependent band-to-band electron-hole 
recombination. The formalism, which is a generaliza- 
tion of a previous spin and charge drift-diffusion theory 
0, applies to both unipolar and bipolar cases, the 
former being a subclass of the latter. Second, we apply 
the formalism to explain the main qualitative features 
of spin accumulation in the extrinsic SHE in the opti- 
cal orientation experiment, for two different boundary 
conditions: (a) uniform generation of electron-hole pairs, 
and (b) edge generation of nonequilibrium electrons. In 
both cases spin accumulation throughout the sample is 
calculated analytically. We find that spin-orbit interac- 
tion modifies the carrier and spin diffusion lengths and 
that the spin accumulation profile depends, qualitatively, 
on the specific boundary conditions, implying that inter- 
pretation of extrinsic SHE requires detailed case-by-case 
considerations for specific experimental and sample ge- 
ometries. 

Consider spin-polarized transport in an inhomoge- 
neous nonmagnetic semiconductor in the presence of elec- 
tric field E. If spin-orbit coupling is present, causing 
skew scattering and side jump, the phenomenological ex- 
pression for the carrier [c — n for electrons and c — p for 



holes) charge current density in the i*'' direction is readily 
obtained by generalizing the Dyakonov-Perel' prediction 

m 

JcX,i = q^icxcxEi ± qDcxdiCx + qXvcXCx<^i]zEj 

+qXdcxeijzdjC\. (1) 

Here the upper (lower) sign is for electrons (holes) and 
A is the spin index; q is the proton charge. The first two 
terms are conventional (longitudinal) carrier drift and 
diffusion, respectively, with /i and D denoting the spin- 
dependent mobility and diffusivity. The third (fourth) 
term represents the effect of skew spin-orbit scattering 
and side jump on drift (diffusion). The effects of the scat- 
tering are in the transverse direction to E and are oppo- 
site for spin up and down carriers. (Holes are treated here 
as spin doublets, which is appropriate in low-dimensional 
structures with heavy and light hole band splitting; oth- 
erwise hole spin does not matter, as we will argue below). 
The corresponding transport parameters are transverse 
mobility v and transverse diffusivity S; they are propor- 
tional to the spin-orbit coupling strength. 

It is more illuminating to introduce the charge, Jc — 
J| + Ji, and spin, = J| — J|, currents. In terms of 
carrier (c = +ci) and spin (sc = Cf — C|) densities, the 
currents are 

Jc,i = q{p-cC + fJ.scSc)Ei ± q{DcdiC + DscdiSc) 

+qeijzEj{vcSc + i^scc) + qeijz{5cdjSc + 5scd;jcl^) 

Jsc,t = 9(MscC + fJ.cSc)Ei ± q{DcdiSc + Dscd^c) 

+qe^j^Ej{vcC + VscSc) + q^tjz{ScdjC + SscdjScjS) 

The transverse carrier charge and spin mobilities are 
given respectively by Vc = {vc-\ + Vci)/'^ and Vsc = 
(/ict~^'cj,)/2, while the transverse carrier charge and spin 
diffusivities are 5c ~ {5c^ +5ci)/2 and 6sc — {5c^ — (5cx)/2. 
The corresponding longitudinal quantities are defined 
similarly. 

Eqs. lO and © succinctly describe the appearance 
of the transverse spin drift and diffusion in the presence 
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of longitudinal charge transport, which is the essence of 
the extrinsic SHE. If the longitudinal current is spin- 
polarized, the above equations describe the anomalous 
Hall effect and the appearance of the transverse Hall 
voltage. 

To further develop the formalism, we need to include 
electron-hole recombination and spin relaxation. In the 
presence of spin-orbit coupling, spin-dependent selection 
rules for band-to-band transitions need to be consid- 
ered. In general, the continuity equation reads: 

±diJcx,t/q = Bi {cxcx ~ cqxCqx) 

+ B2 (cxc^ - coaCqa) + (cA - ci) /2Tit(.4) 

Here c is pin) if c is n^p), cq is the equilibrium carrier 
density, and Tic is the Ti time for spin flipping. 

The spin-preserving recombination rate coefficient Bi , 
as well as the spin-flip coefficient Bo, can be calculated 
by generalizing the unpolarized case |llUl2l| . The valence 
band of zinc-blende semiconductors consists of three sub- 
bands: heavy-hole, light hole, and split-off hole bands. 
We neglect the split-off band as the energy splitting 
A ^ at temperatures T lower than or around room 
temperature. By explicitly taking into account the an- 
gular momentum of the heavy hole and light hole states, 
and using the optical selection rules for the states [lo| . 
we arrive at the following expressions for spin-conserving 
and spin flip recombination constants: 
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(5) 
(6) 



where C depends on the "maximum" electron energy 
^max = Eg + ksT /2 (cg IS the energy gap) and tem- 
perature: 
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C(fta;^ax, T) = j " P'urhu;. (7) 

Here P is the momentum matrix element for optical tran- 
sitions, rij. is the refractive index, TOc (jUh and to/) is the 
band mass of electrons (heavy holes and light holes). 

Eqs. Q and @, together with Poisson's equation form 
a closed set of nonlinear equations whose solution deter- 
mines charge and spin densities and currents in a semi- 
conductor with spin-orbit scattering included. In gen- 
eral, these equations need to be solved numerically for 
specific cases of interest. Our next goal is to introduce 
qualitative features of spin accumulation in two cases 
of experimental interest that allow analytical solutions 
and form a starting point to discuss the concepts and 
issues to be encountered in more complex situations in- 
volving SHE and spin-orbit coupling effects in transport. 
We consider a p-type semiconductor with nondegenerate 



electron (minority) density induced optically. The result- 
ing spin accumulation via extrinsic SHE can be deduced 
in a manner similar to spin orientation experiments. We 
assume that the injected electron density is well below 
the donor density. Further simplification follows from the 
fact that mobilities and diffusivities are spin independent 
in the nondegenerate regime 7J. Finally, we assume un- 
polarized holes since hole spin relaxation in zinc-blende 
semiconductors is extremely fast [l3j. The only carriers 
of interest are then spin-polarized electrons. 

Using the above assumptions, Eqs. 0-0, give the 
drift-diffusion equations for electron spin and carrier den- 
sity: 



{s\/ ■ E + E ■ \/s 



{wp + 1/Tin) S = 0, 



(8) 



(nV ■ E + E ■ Vn) 



(Vs X E)^ 



+w {np - nopo) = 0, (9) 



where w = {Bi -|-i?2)/2 and C = ^n/fJ-n = ^n/Dn charac- 
terizes the spin-orbit coupling strength. The spin quan- 
tization axis is taken to be z. We take E = Ey and 
consider x to be the transverse direction, the slab bound- 
aries being at a; = and a; = a, so that all the quanti- 
ties of interest will have x dependence only. Denoting 
the electron recombination time as Tn = l/{wNa), where 
Na is the acceptor density, and spin relaxation time as 
T, = 1/(t-i + T^^), Eqs. ® and ® become 

Ss/dx"^ + {CqE/kBT)dn/dx ~ s/LI = 0(10) 
d^n/dx^ + {CqE/kBT)ds/dx - {n - no)/Ll = 0(11) 

where the longitudinal spin and charge diffusion lengths 
are defined as Lg = ^DgTs and L„ = V-DnTn, respec- 
tively. In deriving the above equations we have neglected 
terms of order s^, s5n, and Sn"^ relative to 5n — n — uq 
Spin-charge coupling in Eqs. H10|) and (|ll|l is appar- 
ent through the derivatives of spin and charge densities. 

We first solve the transverse spin and carrier diffu- 
sion for the case of a uniformly illuminated slab (bound- 
ary condition BC I), with electron- hole spin- unpolarized 
generation rate G. We assume that carrier recombi- 
nation at the edges is not significant so that it is rea- 
sonable to impose a uniform (here zero) electron trans- 
verse current: Jnx{x) — 0. As for spin current, we take 
Js{x = 0) = Js{x = a) = 0, implying that spin-flip 
scattering at the edges is moderate. The solution to the 
drift-diffusion equation H10|) in the presence of spin-orbit 
scattering is 



s{x)=C 



qELs 
knT 



GTnSec\\{a/2Ls) smh{x/L - a/2L). (12) 



Spin polarization profile is given by a(x) = s{x)/GTn, 
since Gtu ^ no is the average electron density gener- 
ated in steady-state conditions under illumination. This 
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FIG. 1: Calculated spin density s, electron density n, and 
polarization a for BC II. We have used C = 10~* and ap- 
plied electric field strengths 0.125 kV/cm, 0.250 kV/cm, and 
0.500 kV/cm. While a signature of the extrinsic SHE is the 
opposite spin accumulation at the edges of the sample, spin 
accumulation is large also inside. 



simple solution demonstrates the essential physics behind 
spin accumulation in extrinsic SHE: (i) Spin accumula- 
tion increases linearly with E, with possible slight electric 
field modulation due to the dependence of Ls = Ls{E) 
(not discussed here); (ii) The magnitude of spin polar- 
ization is proportional to the strength of the spin-orbit 
scattering as well as to the ratio of the voltage drop over 
min(is, a) and thermal energy; (iii) For a <C is the accu- 
mulation at the edges is linearly proportional to a, while 
for a Ls, the accumulation is independent of a. Fi- 
nally, (iv), while a{x) ^ x for a <ti Lg, spin accumulation 
is significant only within the spin diffusion length from 
the edges when Ls- 

A question now arises: How universal (i.e. indepen- 
dent of boundary conditions) is the qualitative behaviour 
discussed above? To answer this question we introduce 
different boundary conditions (BC II), describing the 
physics of carrier injection at x = and x — a, while 
assuming vanishing spin currents at the edges: n{x = 
0) = ni,n{x = a) = 712; Ja,x{x = 0) = Js,x{x = a) = 0. 
These conditions mimic the case of a p-doped base in a 
pup spin-polarized transistor |l5|. where the electron in- 
jection level can be controlled by the biases to the emitter 
and collector. In the case discussed here the longitudinal 
current would flow perpendicular to the transistor cur- 
rent, which would thus be spin-polarized due to the ex- 
trinsic SHE. The vanishing boundary conditions for Js^x 
reduce to 
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dx 
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(13) 



tions gives the spin and electron densities inside the slab: 
{^qE 



s{x) = 



— cosech( — ) 
Li Li 

cosech( — ) 

L2 L2 



and 

n{x) — riQ 



niCOsh( 



a — x^ 
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Li 
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(14) 
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(72 - 72)cosech(^) 



-(J^_i^)cosech(-^) 



(X — ~ X X 

riisinh(— — ) + n2sinh( — ) 
Li Li 

Qj - — X X 

riisinh(— - — ) + 7i2sinh( — ) 
L2 L2 



(15) 



Here we introduce new transverse spin diffusion lengths, 
Li and L2: 



ir.2 = 7± 



^fF^l/^Lnf 



where 



7 = (1/2) 1/Ll + 1/Ll + (CqE/kBTY 



(16) 



(17) 



Solving Eqs. (|10|l . Hll|) with the above boundary condi- 



There is a critical value of the field which separates 
the regimes of strong and weak spin-charge coupling — 
Eqs. H1U|I and Hill) reduce to the ordinary spin and 
charge diffusion equations when E <ti Es, En , where 
Es — ksT / [C,qLs) and En — ksT / {C,qLn) are the values 
of the critical fields with respect to spin and charge dif- 
fusion. In this case 7 ~ [L~'^ + L~^)/2 and Li^2 — is,n- 
When E > Eg, En, 7 becomes dependent on the elec- 
tric field and it is in this regime the spin-field relation 
deviates from linearity. 

For quantitative understanding we take our semicon- 
ductor to be GaAs at room temperature 0| , with doping 
density Na = 3 x 10"'^^ cm~^, and transverse size of a = 6 
fim (which is much greater than Lg). The spin-orbit scat- 
tering strength is taken to be C = 10^'', reflecting weak 
spin-orbit coupling in GaAs. Finally, the boundary con- 
ditions for electron density are n(0) — 2 x 10^'^/cm'^ and 
n{a) ~ 5x 10^'^/cm'^. Figure Q shows the profiles of 
spin and electron densities, as well as spin polarization 
a = s/n, for several values oi E. In contrast to the purely 
diffusive behavior exhibited by n, spin density along the 
slab is weakly oscillatory. Both spin density and spin 
polarization attain maximum magnitudes at the edges. 
What is interesting is that, unhke in BC I [conclusion 
(iv)], spin accumulation here is significant throughout 
the sample, not only within Ls from the edges. (We 
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FIG. 2: Calculated spin current Js in x (top) and y (bottom) 
directions for BC II and parameters as in Fig. Q 



find that spin polarization is enhanced by a decade when 
considering the extrinsic SHE at 77 K.) 

Spin current is not conserved and flows inside the sam- 
ple, being restricted to zero by our choice of boundary 
conditions. From Eq. Q we obtain the spin current den- 
sity in the x direction: 



< cosech^ 



qVnE 



d ~ — X X 
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and in the y direction: 
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Li 
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Li Li 



— cosech( — ) 
L2 L2 
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L2 L2 



(19) 



Spin current in the x direction flows in the direction of 
the spin gradient, while the y component changes sign 
inside the slab, reflecting the fact that spin-up and spin- 
down electrons are deflected into opposite directions due 
to spin-orbit scattering by impurities. 

Finally, we wish to see at which electric field spin and 
carrier diffusion lengths Li,2 are modified and induce 
nonlinear behavior in a. Figure |21 shows a{E) for differ- 
ent Q. Spin polarization varies linearly with E, except at 



10" 



10' 



10^ 



10- 



E (V/m) 



FIG. 3: Calculated spin polarization a at a; = as a function 



of electric field E for BC II and C = 10" 
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and 1. Only 



in the last (unphysical) case a starts to saturate for large E, 
as the dependence of Li and L2 on E sets in. 



electric fields as large as 10'"' kV/cm and spin-orbit cou- 
plings C « 1, a clearly unphysical case considered here 
only to illustrate the scope of linear behavior. 

In conclusion, we have presented a drift-diffusion for- 
malism which takes into account spin-orbit scattering 
and spin-dependent carrier recombination. We have 
calculated spin accumulation in the extrinsic spin Hall 
regime and introduced spin-orbit dependent carrier and 
spin diffusion lengths. We have found that spin accu- 
mulation is strongly influenced by boundary conditions 
and thus by specific spintronic device design. Our the- 
ory is applicable to any complex device setting in which 
extrinsic spin Hall effect is expected to play a role. We 
expect similar strong dependence of the "intrinsic" SHE 
on boundary conditions too, making it difficult to distin- 
guish intrinsic and extrinsic SHE in general. 

This work was supported by the US ONR, NSF, and 
National Research Council (I.Z.). 
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